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Abstract 

The solar cell structure of n -type poly-silicon/5-nm-diameter silicon nanocrystals embedded in an amorphous silicon 
oxycarbide matrix (30 layers)/p-type hydrogenated amorphous silicon/AI electrode was fabricated on a quartz substrate. 
An open-circuit voltage and a fill factor of 518 mV and 0.51 in the solar cell were obtained, respectively. The absorption 
edge of the solar cell was 1.49 eV, which corresponds to the optical bandgap of the silicon nanocrystal materials, 
suggesting that it is possible to fabricate the solar cells with silicon nanocrystal materials, whose bandgaps are wider 
than that of crystalline silicon. 
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Background 

Over the past few years, many researchers have shown 
an interest in silicon nanostructures, such as silicon 
nanocrystals [1-4] and silicon nanowires [5-8] for solar 
cell applications. Since a silicon nanocrystal embedded 
in a barrier material can make carriers confined three- 
dimensionally, the absorption edge can be tuned in a 
wide range of photon energies due to the quantum size 
effect. Thus, it is possible to apply silicon nanocrystal 
materials or silicon quantum dot (Si-QD) materials to all 
silicon tandem solar cells [9], which have the possibility 
to overcome the Shockley-Queisser limit [10]. Moreover, 
it has been found that the weak absorption in bulk Si is 
significantly enhanced in Si nanocrystals, especially in the 
small dot size, due to the quantum confinement-induced 
mixing of T-character into the X-like conduction band 
states [11]. Therefore, Si-QD materials are one of the 
promising materials for the third-generation solar cells. 
Size-controlled Si-QDs have been prepared in an amorph- 
ous silicon oxide (a-Si0 2 ) [12], nitride (a-Si 3 N 4 ) [13], carbide 
(a-SiC) [14-17], or hybrid matrix [18,19], which is called as 
silicon quantum dot superlattice structure (Si-QDSL). In 
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the case of solar cells, generated carriers have to be 
transported to each doping layer. Since the barrier height 
of an a-SiC matrix is relatively lower than that of an a- 
Si 3 N 4 or a-Si0 2 matrix, the Si-QDSL using an a-SiC 
matrix has an advantage in carrier transport. Therefore, 
the development of the Si-QDSL solar cells using an a-SiC 
matrix is of considerable importance. There are a few 
researches fabricating Si-QDSL solar cells. Perez- Wurfl 
et al. reported that Si-QDSL solar cells with Si0 2 matrix 
showed an open-circuit voltage ( V oc ) of 492 mV. However, 
the clear evidence of the quantum size effect has not been 
reported from Si-QDSL solar cells [20]. In our previous 
work, Si-QDSLs with a-SiC matrix have been prepared by 
plasma-enhanced chemical vapor deposition (PECVD). 
The defect density in a Si-QDSL has been successfully re- 
duced by hydrogen plasma treatment (HPT), and the shift 
of photoluminescence spectra from the Si-QDSLs has been 
confirmed by varying the diameter of the Si-QDs [2]. More- 
over, it has been revealed that the oxygen-incorporation 
into the a-SiC matrix can suppress the formation of the 
leakage paths [21]. An V oc of 518 mV has been obtained in 
a Si-QDSL solar cell with an amorphous silicon oxycarbide 
(a-Sii _ x _ yC x O y ) matrix [1]. 

In this paper, we report the effect of oxygen addition 
on the formation of Si-QDs in a-Sii _ x _ y C x O y . Optical ab- 
sorption coefficients of the Si-QDSL were also investigated. 
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Si-QDSL solar cells were fabricated using the optimum 
oxygen concentration. In addition, the numerical analysis 
using the Bohm quantum potential (BQP) method was 
performed to simulate the electrical characteristics of fabri- 
cated solar cells. 



Methods 

Experimental method 

The a-Sii _ x _ y C x O y matrix was deposited on a quartz sub- 
strate to investigate the fundamental optical properties such 
as Raman scattering spectrum, transmittance, and re- 
flectance. The fabrication method is referred as follows. 
A 40-period-multilayer with silicon-rich hydrogenated 
amorphous silicon oxycarbide layers and hydrogenated 
amorphous silicon oxycarbide barrier layers was prepared 
on a quartz substrate by very high frequency PECVD 
method (VHF- PECVD). The source gases were silane 
(SiH 4 ), monomethylsilane (MMS), hydrogen (H 2 ), and car- 
bon dioxide (C0 2 ). The flow rates of SiH 4 , MMS, and H 2 + 
C0 2 ; deposition pressure; substrate temperature; frequency; 
and plasma power were fixed at 3.3 , 1.3, and 47.4 seem; 
20 Pa; 60 MHz; 193 °C; and 13 mW/cm 2 , respectively. The 
flow rate of C0 2 was varied from 0 to 3.7 seem. The mass 
flow controllers for SiH 4 and C0 2 were calibrated by N 2 . A 
H 2 -calibrated mass flow controller was used for MMS. 
During the deposition of a-Si 2 _ x _ y C x O y barrier layers, the 
flow of SiH 4 gas was stopped. Subsequently, the samples 
were annealed at 900 °C for 30 min under a forming gas 
atmosphere to form Si-QDs in an a-Si 2 _ x _ y C x O y 
matrix. The target size of Si-QDs and barrier width were 
5 and 2 nm, respectively. The concentrations of Si, C, and 
O in the barrier layer were measured by X-ray photoelec- 
tron spectroscopy (XPS). The crystallinity of Si-QDs was 
investigated by Raman scattering spectroscopy. The ab- 
sorption coefficient of a Si-QDSL was estimated by the 
transmittance and the reflectance of a sample. The samples 
with uniform thickness were selected for the measure- 
ments, and one measurement was carried out for each 
measurement method and for each sample. 

The solar cells using Si-QDSL as an absorber layer were 
also fabricated. The schematic of the solar cell structure is 
shown in Figure 1. The fabrication process is referred as 
follows. A phosphorus-doped hydrogenated amorphous 
silicon thin film was deposited on a quartz substrate by 
PECVD. The film was annealed at 900°C for 30 min under 
a forming gas, resulting in a polycrystalline silicon (n-type 
poly-Si) thin film. On the poly-Si layer, a 30-period super- 
lattice was deposited by VHF-PECVD. The same depos- 
ition conditions as referred to above were used except the 
flow rates of C0 2 and H 2 , which were fixed at 0.4 and 
47.0 seem, respectively. Consequently, the sample was 
annealed at 900°C for 30 min to form Si-QDs. The sample 
was exposed to hydrogen plasma to reduce dangling bond 
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Figure 1 The schematic of the fabricated Si-QDSL solar 
cell structure. 



defects in the post-annealed Si-QDSL. After the treat- 
ment, a boron-doped hydrogenated amorphous silicon (p- 
type a-Si:H) was deposited by PECVD. An aluminum (Al) 
electrode was deposited by the evaporator on the sample. 
The electrode area of a solar cell was 0.00785 cm 2 . The 
cross-sectional structure of a solar cell was observed by 
transmission electron microscopy (TEM). The solar 
cells were characterized by dark I-V characteristics and 
light I-V characteristics under the illumination at 
AM1.5G, 100 mW/cm 2 , and 25°C. 



Numerical method 

The numerical calculations of the Si-QDSL solar cells 
were performed using a two-dimensional device simulator, 
Atlas ver. 5.18.3.R (Silvaco, Inc., Santa Clara, CA, USA). 
The device structure used for numerical calculations is 
shown in Figure 2. Quartz substrate/n-type poly-Si/30- 
period Si-QDSL (Si-QDs embedded in a-Sii _ x _ y C x O y )/p- 
type hydrogenated amorphous silicon (p-type a-Si:H)/Al 
electrode structure was assumed in this simulation. The 
diameter of Si-QDs and the gap between any two Si-QDs 
were fixed at 5 and 2 nm, respectively. The BQP method 
[22-26] was adopted to describe the quantum confine- 
ment effect and the quantum tunnel effect in the Si-QDSL 
layer. The electrical transport in the Si-QDSL was 
described by drift- diffusion equations and current con- 
tinuity equations for electrons and holes. In the theory, 
the transport of carriers is influenced by the total potential 
of the potential characterizing the system and quantum 
potential. The definition of the effective quantum poten- 
tial Q e ff is derived from a weighted average of the BQPs 
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Figure 2 The structure of the Si-QDSL solar cell for 
numerical calculations. 



seen by all single-particle wavefunctions, which can be 
expressed as 



Qeff, n — ~Y n 



h 2 v 2 y^ 

8jT 2 m n ^Jn 



and 



<?eff, p — ~Y J 



h 2 V 2 y^ 
p 8n 2 rn p ^/p 



(i) 



(2) 



where Q e ff,« and Q e ff, p are the effective BQPs for the con- 
duction band and the valence band, respectively, h is 
Planck's constant, n and p are the electron and hole con- 
centrations, respectively. y n and y p are adjustable param- 
eters for quantum confinement. In general, a three- 
dimensional quantum system cannot be described on a 
two-dimensional device simulator due to the difference 
of the quantum confinement effects between two- and 

Table 1 Parameters of each layer for calculations 

Parameters n-type Si-QD a-Sh _ x _ y C x O y p-type 





poly-Si 






a-Si 


Energy gap (eV) 


1.13 


1.13 


2.5 


1.7 


Electron affinity (eV) 


4.17 


4.17 


3.5 


4.0 


Carrier lifetime (s) 


1 X10" 15 


1 x 1CT 10 


1 x 10~ 10 


1 x 10~ 6 


Electron mobility 

(cm 2 /Vs) 


1 


1 


1 


1 


Hole mobility 

(cm 2 /Vs) 


0.1 


0.1 


0.1 


0.1 


Donor concentration 

(cm" 3 ) 


1 X10 19 


0 or 

1 x10 17 






Accepter concentration 

(cm" 3 ) 








1X10 19 



Table 2 Concentrations of Si, C, and O in a-Si n . 
films with several C0 2 /MMS flow rate ratios 



C0 2 /MMS 


Si (at.%) 


C (at.%) 


O (at.%) 


0 


44.6 


37.9 


17.5 


0.3 


40.3 


34.6 


25.1 


1.5 


34.2 


33.2 


32.6 


3.0 


31.9 


28.3 


39.8 



three-dimensional systems. To take three-dimensional 
quantum effect into two-dimensional simulation, we ad- 
justed the Yn an d Yp parameters in the BQP model. The 
parameter values were determined to satisfy that the 
bandgap calculated from the BQP method is equal to 
the bandgap derived from three-dimensional Schrodin- 
ger equations. The y n and the y p for the Si-QDSL with 
5-nm-diameter Si-QDs and 2-nm-thick a-Sii _ x _ y C x O y 
barrier layers were 4.0. In this simulation, the radiative 
recombinations, Shockley-Read-Hall (SRH) recombin- 
ation [27-29] and Auger recombination, were taken into 
account. Auger coefficients and effective masses of bulk Si 
were adapted for all layers. The other parameters are 
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Figure 3 The Raman spectra of the Si-QDSLs with several C0 2 / 
MMS flow rate ratios, (a) C0 2 MMS = 0. (b) C0 2 MMS = 0.3. 
(c) C0 2 MMS = 1 .5. (d) C0 2 MMS = 3. 
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shown in Table 1. The bandgaps in the table do not affect 
optical absorption but carrier transport phenomenon. To 
take into account the phosphorus diffusion into the Si- 
QDSL layer, a calculation with the donor concentration in 
Si-QDs of 1 x 10 17 cnT 3 was also performed. The light I-V 
characteristics were calculated, assuming solar illumin- 
ation of AM1.5G at 100 mW/cm 2 . Additionally, the 
quantum efficiencies were calculated without bias light 
and bias voltage. An incident light was put into the solar 
cells from the quartz substrate side normally. The light in- 
tensity and the photogeneration rate were calculated 
based on the ray tracing method, where the Si-QDSL was 
regarded as an optically homogeneous material, and the 
optical parameters from the spectroscopic ellipsometry 
measurement of the Si-QDSL were used. 




Photon energy (eV) 

Figure 4 The absorption coefficients of the Si-QDSLs with several 
C0 2 /MMS flow rate ratios. 

V J 



Results and discussion 

Optical properties of Si-QDSLs 

The concentrations of Si, C, and O in a-Sii _ x _ y C x O y 
thin films were measured by the relative sensitivity fac- 
tor (RSF) method. The concentrations of Si, C, and O 
for each C0 2 /MMS flow rate ratio were shown in 
Table 2. The oxygen concentration and the deposition 
rate of the films depend on the C0 2 /MMS flow rate 
ratio. The oxygen concentrations of the films prepared 
without C0 2 gas and with the C0 2 /MMS flow rate 
ratios of 0.3, 1.5, and 3.0 were 17.5, 25.1, 32.6, and 
39.8 at.%, respectively. Oxygen was observed even in the 
as-deposited film prepared without flowing C0 2 gas. 
This unintentionally incorporated oxygen is thought to 
be originating from the deposition atmosphere. The de- 
position rate is proportional to the oxygen concentration 
in the film, suggesting that the volume of the thin film 
increases with the oxygen incorporation. 

The crystallization of Si-QDs was investigated by Raman 
scattering spectroscopy. The Raman spectra of the Si- 



QDSLs with the C0 2 /MMS flow rate ratios of 0, 0.3, 1.5, 
and 3.0 are shown in Figure 3. A Raman spectrum was 
separated into three Gaussian curves. The peaks at ap- 
proximately 430 and 490 cm" 1 are originating from the 
LO mode and TO mode of a-Si phase, respectively [30]. 
These Gaussian curves are indicated by blue dashed lines. 
The peak at approximately 510 cm" 1 is originating from 
Si-QDs. The Gaussian curve is indicated by green dashed 
line. As the C0 2 /MMS flow rate ratio increases, the inten- 
sity of the peak from Si-QDs becomes weaker compared 
with the peak from a-Si phase. This indicates that the 
crystallization of Si-QDs in the silicon-rich layers is pre- 
vented by the oxygen-incorporation, and the crystallization 
temperature of nanocrystalline silicon phase becomes 
higher [31]. 

The absorption coefficient was estimated from the mea- 
surements of transmittance and reflectance. The absorp- 
tion coefficients of the Si-QDSLs with the C0 2 /MMS flow 
rate ratios of 0, 0.3, 1.5, and 3.0 are shown in Figure 4. For 
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Figure 5 The cross-sectional TEM images of the fabricated solar cell structure, (a) The whole region image with the schematic of the 
structure and the thicknesses of each layer, (b) The magnified image of the Si-QDSL layer in the solar cell. 
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both Si-QDSLs with the C0 2 /MMS flow rate ratios of 0 
and 03, the absorption enhancement was observed below 
the photon energy of 2.0 eV. Moreover, the absorption 
enhancement becomes weaker as the C0 2 /MMS flow rate 
ratio increases. This tendency corresponds to that of the 
intensity of the peak originating from Si-QDs in the 
Raman scattering spectrum. Therefore, one can conclude 
that the absorption enhancement is due to the increment 
of the nanocrystalline silicon phase. Moreover, the absorp- 
tion edge was estimated by the Tauc model [32]. The 
absorption edges of the Si-QDSLs with the C0 2 /MMS 
flow rate ratios of 0 and 0.3 were estimated at 1.48 and 
1.56 eV, respectively. These values are similar to the op- 
tical gap of 5-nm-diameter Si-QDs in an a-SiC matrix 
measured by photoluminescence spectrum [2]. On the 
other hand, the absorption edges of the Si-QDSLs with 
the C0 2 /MMS flow rate ratios of 1.5 and 3.0 were esti- 
mated at approximately 1.70 eV, which corresponds to the 
optical gap of a-Si. 

These results indicate that the C0 2 /MMS flow rate 
ratio should be below approximately 0.3 to form 
Si-QDs in the silicon-rich layers. According to the 
[22], the CO2/MMS flow rate ratio should be higher 
than 0.3 to suppress the crystallization of a-SiC phase 
in the a-Sii _ x _ y C x O y barrier layers and the increment 
of the dark conductivity for the annealing temperature 
of 900°C. Although there is a trade-off between the 
promotion of the crystallization of Si-QDs and the 
suppression of the crystallization of a-SiC phase, the 
C0 2 /MMS flow rate ratio of approximately 0.3 or 
the oxygen concentration of approximately 25 at.% is 
one of the optimal conditions. Therefore, the C0 2 / 
MMS flow rate ratio of 0.3 is adopted for the solar cell 
fabrication in this study. 



light I-V (CO 2 /MMS=0.3) 
Dark I-V (CO 2 /MMS=0.3) 
light I-V (CO 2 /MMS=0) 
Dark I-V (CO 2 /MMS=0) 




0 0.2 0.4 

Voltage (V) 

Figure 6 The I-V characteristics of the fabricated Si-QDSL solar 
cell [1,3]. 



Table 3 Solar cell parameters of the fabricated Si-QDSL 
solar cells and the calculated by BQP method 



Parameters 


Experimental 


Calculated 






Doped Si-QDSL Non-doped Si-QDSL 


V oc (mV) 


518 


520 505 


J sc (mA/cm 2 ) 


0.34 


3.98 4.96 


FF 


0.51 


0.61 0.69 



I-V characteristics of the fabricated solar cells 

The cross-sectional TEM images of the fabricated solar 
cell are shown in Figure 5. Figure 5a shows the image of 
the whole region of the solar cell. Figure 5b shows the 
magnified image of the Si-QDSL layer in the solar cell. 
The thicknesses of the n-type poly-Si layer, the Si-QDSL 
layer, and p-type a-Si:H layer were approximately 530, 
143, and 46 nm, respectively. The black region below 
the n-type poly-Si layer is a quartz substrate. The tex- 
tured quartz substrate is used to prevent from peeling 
off the films during the thermal annealing. In Figure 5b, 
the yellow lines and orange circles indicate the interface 
between an a-Sii _ x _ y C x O y barrier layer and a Si-QD 
layer, and Si-QDs, respectively. This magnified image 
revealed that a Si-QDSL layer including average 5-nm- 
diameter Si-QDs was successfully prepared. 

Figure 6 shows the dark I-V characteristics and the 
light I-V characteristics of the solar cells with the C0 2 / 
MMS flow rate ratio of 0 and 0.3 [1,3]. The diode prop- 
erties were confirmed from the dark I-V characteristics. 
The characteristics were evaluated by one-diode model: 



: I 0 exp 



nkT 



V + RsI 



(3) 



where I 0 , n, R s , and R sh represent reverse saturation 
current density, diode factor, series resistance, and shunt 
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Figure 7 Internal quantum efficiencies of fabricated solar cell 
and of that calculated by the BQP method. 
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resistance, respectively. According to the fitting of the 
dark I-V characteristics of the oxygen-introduced Si- 
QDSL solar cell, the reverse saturation current density, 
the diode factor, the series resistance, and the shunt 
resistance were estimated at 9.9 x 1(T 6 mA/cm 2 , 2.0, 
2.3 x KT 1 a cm 2 , and 2.1 x 10 4 H cm 2 , respectively. The 
solar cell parameters of the light I-V characteristics 
under AM1.5G illumination are summarized in Table 3. 
An V oc of 518 mV was achieved. Compared with the V oc 
of 165 mV with non-oxygen-introduced Si-QDSL solar 
cells, the characteristics were drastically improved. The 
possible reasons for this improvement are due to the 
passivation effect of Si-O phase on silicon quantum dots 
[33], and the reduction of the leakage current by the 
introduction of oxygen [21]. Figure 7 shows the internal 
quantum efficiency of the solar cell. The red line corre- 
sponds to the experimental internal quantum efficiency. 
The quantum efficiency decays to zero at approximately 
800 nm, suggesting that the contribution is originating 
not from the n-type poly-Si but from the Si-QDSL ab- 
sorber layer. 

Calculations of I-V characteristics and quantum 
efficiencies 

The light I-V characteristics and the internal quantum effi- 
ciency of the Si-QDSL solar cells with a doped Si-QDSL 
layer and a non-doped Si-QDSL layer were simulated 
under AM1.5G illumination using the BQP method. The 
calculated solar cell parameters are shown in Table 3. 
Also, the calculated quantum efficiencies are shown in 
Figure 7. The simulated quantum efficiencies are multi- 
plied by 0.12 for comparison with the experimental one. 
The calculated short-circuit current densities (J sc ) and 
quantum efficiencies are much higher than those of the 
experimental results. There are two possible reasons. 

The first reason is due to the difference of the doping 
concentration in a Si-QDSL layer. In an actual solar cell, 
the phosphorus concentration in the Si-QDSL absorber 
layer is more than 1 x 10 19 cm" 3 due to the high- 
temperature annealing process [34]. From the simula- 
tions, the and the quantum efficiency in the whole 
wavelength region becomes lower if the phosphorus 
concentration in the Si-QDSL layer increases. The phos- 
phorus in the Si-QDSL layer degrades the / sc due to the 
reduction of the electrical field in the Si-QDSL layer. 
Unfortunately, simulations were not possible when the 
dopant concentration in the Si-QDSL was higher than 
1 x 10 17 cm" 3 due to the convergence problem of the 
BQP calculations. It is expected that / sc will decrease 
more if the dopant concentration becomes higher. We 
previously reported that the quantum efficiency in the 
whole wavelength region decreases as the dopant con- 
centration in the Si-QDSL increases from experiments 
and the simulations using classical model [35], which is 



similar to the results of the BQP method. The second 
reason is due to the optical losses in the n-type poly-Si 
layer. In this calculation, the surface roughness of the 
textured quartz substrate was not taken into account. 
The effective optical path length in the n-type layer of 
the simulated structure should be shorter than that of 
the actual solar cell structure. As a result, the simulated 
quantum efficiency in the short-wavelength region is 
higher than that of the experimental because of the low 
optical absorption loss in the n-type poly-Si layer. 

Even though the / sc mismatch, the absorption edge can 
be estimated from the simulated quantum efficiency. The 
calculated quantum efficiencies at the long-wavelength re- 
gion are in agreement with those of the experimental one. 
This suggests that the absorption edge of the solar cell can 
be theoretically reproduced using this simulation. More- 
over, the absorption edge was estimated to be 1.49 eV, 
which is quite similar to the absorption edge of the Si- 
QDSL estimated from the optical measurements. This in- 
dicates that the photogeneration in the Si-QDSL solar cell 
is thought to be the contribution from Si-QDs, and it is 
possible to fabricate the solar cells with silicon nanocrystal 
materials, whose bandgaps are wider than that of a crystal- 
line silicon. 

Conclusions 

The fundamental optical properties of Si-QDSLs were in- 
vestigated, and the solar cell structure using the Si-QDSL 
as an absorber layer was fabricated and characterized. 
From the measurements of the Raman spectra and the ab- 
sorption coefficients of Si-QDSLs, it was revealed that the 
absorption coefficient is enhanced by the crystallization of 
the Si-QDs, and the crystallinity of Si-QDs is affected by 
the oxygen concentration in the superlattice. In addition, 
the solar cell characteristics were simulated by the BQP 
method. The absorption edge of the simulated Si-QDSL 
solar cell was in agreement with that of the fabricated one. 
Moreover, the absorption edge of the Si-QDSL solar cell 
was 1.49 eV, which is similar to the absorption edge esti- 
mated from the optical measurements. These results sug- 
gest that it is possible to fabricate the solar cells with 
silicon nanocrystal materials, whose bandgaps are wider 
than that of a crystalline silicon. 
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